The Feynman propagator for spin foam quantum gravity 
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We link the notion causality with the orientation of the 2-complex on which spin foam models 
are based. We show that all current spin foam models are orientation-independent, pointing out 
the mathematical structure behind this independence. Using the technology of evolution kernels 
for quantum fields/particles on Lie groups/homogeneous spaces, we construct a generalised version 
of spin foam models, introducing an extra proper time variable and prove that different ranges of 
integration for this variable lead to different classes of spin foam models: the usual ones, interpreted 
as the quantum gravity analogue of the Hadamard function of QFT or as a covariant definition of 
the inner product between quantum gravity states; and a new class of causal models, corresponding 
to the quantum gravity analogue of the Feynman propagator in QFT, non-trivial function of the 
orientation data, and implying a notion of "timeless ordering" . 



I. INTRODUCTION 

Spin foam models Q, are a new more rigorous im- 
plementation of the path integral or sum-over-histories 
approach to quantum gravity, in which spacetime is re- 
placed by a combinatorial 2-complex, and the histories 
of the gravitational field are characterized by data taken 
from the representation theory of the local gauge group of 
gravity, i.e. the Lorentz group, and no familiar notions of 
metric on differentiable manifolds are used. The main ob- 
ject that such a sum-over-histories formulation defines is 
the partition function of the theory, and from this, allow- 
ing the underlying spacetime 2-complex to have bound- 
aries, quantum amplitudes functions of the boundary 
data, that, when a canonical interpretation is available, 
should be interpreted as transition amplitudes between 
quantum gravity states. The partition function in these 
models, for a given 2-complex a, has the general struc- 
ture: 

= E I[Mpf)I[Mpf\e)l[Mpf\v), (i) 

{p/} / e v 

where the amplitudes for faces /, edges e and vertices 
v of the 2-complex are all functions of the representa- 
tions pf of the Lorentz group associated to the faces of 
the complex. This work deals, in the spin foam setting, 
with two related issues common, but not restricted, to 
any sum-over-histories formulation of quantum gravity: 
1) given a definition of transition amplitudes for quan- 
tum gravity states, which kind of transition amplitude 
does it provide? 2) is it possible to identify a notion of 
causality and to implement some causality conditions ex- 
plicitly in quantum gravity, in spite of the necessary (be- 
cause of diffcomorphism invariance) absence of any time 
parameter in the quantities one wants to compute? The 
two issues are related because it is known that in the 
case of free relativistic particle and quantum field the- 
ory the main difference between the different transition 
amplitudes one can construct, in particular the Feyn- 



man propagator and the Hadamard function, is exactly 
the way they encode (or fail to encode) causality restric- 
tions, i.e. the fact that one of their arguments is in the 
causal future of the other, or in other words the order 
between these arguments. This order-dependence can be 
imposed in each of the histories summed over in a clear 
way using a proper time formulation |J 5 : starting from 
the same proper time dependent expression g{x\, X2, T), 
the Hadamard function is obtained by integrating over 
both positive and negative values of the proper time; the 
resulting amplitude is then a-causal, and does not reg- 
ister any ordering between its arguments; the Feynman 
propagator is instead obtained restricting this integra- 
tion to positive proper times only, and this ordering is 
precisely what makes it causal. What about quantum 
gravity? Here no external time coordinate is available, 
so no definition of the Feynman propagator in terms of 
time-ordered products is possible; however, as we said, 
the crucial property that distinguishes the two Green 
functions of the relativistic particle is the symmetry un- 
der a change in the orientation of the history connecting 
the initial and final states and, when a proper time for- 
mulation is available, this is what the signs of the proper 
time amounts to. Therefore we see a way out of the men- 
tioned difficulty in the quantum gravity case. Indeed it 
has been proposed Q that a similar situation as in the 
particle case arises in quantum gravity when formulated 
in the path integral formalism; the distinction between 
the different transition amplitudes is again obtained by 
choosing different range of integration in the proper time 
formulation of the theory; using a canonical decomposi- 
tion of the metric variables, the role of proper time is 
played by the lapse function; again an unrestricted range 
of integration leads to the (quantum gravity analogue of 
the) Hadamard function, while a restriction to positive 
lapses leads to the (analogue of the) Feynman propagator 
The problem we address in this work is how to con- 
struct the analogue quantities in the spin foam context, 
i.e. in a purely combinatorial, algebraic and group the- 
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oretic way, in absence of any smooth manifold structure 
and any metric field. 



II. CAUSALITY AS ORIENTATION 

In spin foam models spacetime is replaced by a combi- 
natorial 2-complex, i.e. by a collection of vertices, links 
and faces, with extra data assigning more geometric in- 
formation to it. What can be the analogue of causal rela- 
tions in such a context? Consider just the first layer of the 
spin foam 2-complex, i.e. only vertices and links connect- 
ing them, therefore forming a graph. If we add to it ori- 
entation data, i.e. arrows on the links, we obtain an ori- 
ented (or directed) graph, a set of oriented links connect- 
ing a set of vertices. Now the vertices can be interpreted 
as a set of fundamental spacetime events and the oriented 
links are then the causal relations between them. We can 
assign an orientation variables a e i v = ±1 to each link, 
with respect to each vertex v it connects, and [i = ±1 
to each vertex. The spacetime interpretation of these 
variables is that of indicating whether the vertex is a fu- 
ture pointing or past pointing contribution to the overall 
spacetime diagram, and whether the link is ingoing or 
outgoing with respect to each vertex. A consistency con- 
dition for the assignment of orientation data to the graph 
is that when a link e connects two vertices it has the 
opposite orientation in the two: a e \ Vl p Vl — —a e \ V2 /j IV2 . 
We also assign an extra orientation variable to each face 
6f = ±1. This structure is basically that of a poset or 
causal set Q , but the set of vertices-events endowed with 
the ordering relation now fails to satisfy in general any 
reflexivity condition, in other words, our causal relations 
allow for closed timelike loops. Also, while the combina- 
torial structure above is general, its causal interpretation 
makes sense only in a Lorentzian context; therefore the 
issue we will be confronting in the following is the gen- 
eral one of constructing spin foam models that take into 
account appropriately the orientation of the underlying 
2-complex, i.e. of orientation- dependent transition am- 
plitudes for quantum gravity. In a Lorentzian context 
these will have the interpretation of causal amplitudes 
or of quantum gravity analogues of the Feynman prop- 
agator. The 2-complexes used in spin foam models are 
not generic: they are topologically dual to simplicial n- 
dimcnsional complexes: to each vertex corresponds a n- 
simplex, to each link a (n-l)-dimensional simplex, to each 
face an (n-2)-dimensional simplex. The orientation data 
we assigned have then a clear geometric interpretation 
in this simplicial picture: the fi v variable for a vertex 
takes the values ±1 according to whether the n-simplex 
dual to it is isomorphic to a n-simplex in Minkowski (Eu- 
clidean) space or the isomorphism holds for the opposite 
orientation; the variable a e \ v = ±1 indicate whether the 
normal to the (n-l)-simplex dual to the link e is ingoing 
or outgoing with respect to the n-simplex dual to v, and 
the variables e/ also characterize the orientation of the 
(n-2)-simplex dual to the face /. A second consistency 



condition on the values of the orientation variables, that 
basically follows from Stokes's theorem, is then 0, @: 
Vw ey|. u = a ei |. u a e2 | 1 ,/i„, where e\ and ei label the two 
links that belong to the boundary of the face / and touch 
the vertex v. This orientation structure is the quantum 
seed of the classical causal structure, in a Lorentzian con- 
text. 



III. ORIENTATION-INDEPENDENCE OF 
CURRENT SPIN FOAM MODELS 

We have now to see how the orientation data intro- 
duced in the previous section enter in current spin foam 
models. Our analysis holds for the Ponzano-Regge [6j 
and all Barrett-Crane-like models Q, based on ordi- 
nary Lie groups and homogeneous spaces, with no ex- 
plicit quantum group structure. Physically these can 
be interpreted as models of quantum gravity without 
cosmological constant. We use a 'first order' formula- 
tion of spin foam models in terms of both group 
variables (or variables with values in an homogeneous 
space) and representation variables. Let us anticipate 
the result before justifying it: all current spin foam mod- 
els do not depend, in their amplitudes, on the orien- 
tation of the underlying 2-complex. The way this is 
achieved is quite simple: in the expression for the am- 
plitudes for spin foams the terms that can be understood 
as contributions from opposite orientations are summed 
simmetrically thus erasing the dependence on the ori- 
entation data. We deal here explicitely only with the 
n-dimcnsional models for n > 4 based on the homo- 
geneous spaces SO(n — 1,1)/ SO(n — 1) ~ and 
SO(n)/SO(n — 1) ~ 5 n_1 for brevity, but the analysis 
can be extended with the same results for all the other 
spin foam models The Barrett-Crane-like models in 
4 and higher dimensions [lj take the form: 

z = (nf d Pfj (nn J H ^ n^fipf) 

e v 

where of course the precise combinatorics varies accord- 
ing to the dimension, but in any case: pf are the unitary 
irreps of the local gauge group of gravity (SO(n — 1,1) 
or SO(n)), and these (class I) unitary representations are 
labelled by either a half-integer j in the Ricmannian case, 
or by a positive real parameter p in the Lorentzian, H e 
is the homogeneous space to which the vectors x e , inter- 
preted as normals to (n-l)-simplices in each n-simplex, 
belong. The vertex amplitudes factorize in terms of face 
contributions, and it is at this level that the information 
about the orientation of the 2-complex is erased: 

A v (x eev ,p fev ) = Y[A fev {9f,p f ) = 
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where the amplitudes depend on the x e vectors only 
through the invariant distances 6j — cosh(/i) (a; e i e / ■ 
a; e 2e/), and in the last step we have traded the orienta- 
tion data on the faces for those on the vertices by us- 
ing the second consistency relation, indicating as df the 
oriented angle a e \a e iQf. This orientation independent 
structure is pretty general (independent on signature or 
dimension), and such are the mathematics behind this 
and the orientation dependent functions W. In fact the 
amplitudes Af^ v are constructed out of representation 
functions of the relevant Lie group, Af^ v = D^^Of), i.e. 
they are given by zonal spherical functions. In turn all 
these representation functions can be uniquely expressed 
in terms of sums of so-called "representation functions 
of the 2nd kind" E?{ (g) ^(ji solutions of the same set 
of equations with different boundary conditions , 
and it is precisely these representation functions of the 
2nd kind that originate the orientation-dependent am- 
plitudes Wf£ V HI- The relation between representation 
functions of the 1st and 2nd kind is best understood in 
terms of their analytic expressions in terms of associated 
Legendre functions of the 1st and 2nd kind 0, ^| : 



2T^r(^i| 3- 

Afev = D p (8 f ) = p a •■ 

(sin(/i)0)M 
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(sm{h)e)^ r (a + 2=1) r (-a - a=2) cos™ 
"QffliW - qT{0)] = E p oo (0) + E m p {9) -- 
Wf= +1 (^,p / )+Wf=- 1 (^ / , P/ ) 



where the a = 2j + (n — 3)/2 in the n-dimensional Ric- 
mannian case, and a — ip — 1/2 in the Lorentzian case. A 
similar relation can be obtained 9] for the 3-dimensional 
Ponzano-Regge case, and for the other Lorentzian mod- 
els. This orientation independence leads to interpret- 
ing the current spin foam models as a-causal transition 
amplitudes, i.e. as the quantum gravity analogue of 
the Hadamard function; when a canonical formulation 
is available, they can be equivalently thought of as defin- 
ing the physical inner product between quantum gravity 
states, invariant under spacetimc diffcomorphisms, or as 
a covariant definition of the matrix elements of the pro- 
jector operator onto physical states Given the uni- 
versal structure outlined above, we are lead to look for 
a universal way of modifying current spin foam models 
and to a new, again universal, definition of orientation- 
dependent or causal spin foam models [j|; these would 
correspond to causal transition amplitudes and to a quan- 
tum gravity analogue of the Feynman propagator of field 
theory. The brute force procedure is obvious and was 
performed in the 4-dimensional Lorentzian case in 0: 
simply restrict the spin foam amplitudes switching from 
D functions to E functions corresponding to a consistent 



choice of values for the orientation data, dropping the 
sum over p v or ey. Hovewer, we will shortly see that a 
much more natural and elegant construction exists. Also, 
the expressions |21 are reminescent of the decomposition of 
Hadamard functions into Wightman functions. The new 
construction we are proposing shows that this analogy is 
indeed exact. 



IV. SPIN NETWORK FEYNMANOLOGY AND 
PARTICLES ON LIE GROUPS/HOMOGENEOUS 
SPACES 

The amplitudes assigned to spin foam faces, edges 
and vertices in the Barrett-Crane-type models is given 
by the evaluation of simple spin networks and was de- 
scribed in [13j in analogy to the evaluation of Feynman 
diagrams: assign a variable valued in the relevant ho- 
mogeneous space to each vertex of the given spin net- 
work, assign a zonal spherical function D^(8f) to each 
line and sum over all the possible values of the vari- 
ables on the vertices to get the final amplitude. In 
this prescription the zonal spherical function is treated 
as a kind of propagator, and this is what it really is, 
as it turns out. Consider a scalar field <p(g) with mass 
to living on the Lie group/homogeneous space G, with 
each point on it labelled by g; consider its free evolution 
parametrised by a proper time coordinate s; the equa- 
tion of motion in proper time is: (id s + A)(f>(g, s) — 
with A being the Laplace-Beltrami operator on G. 
The dynamics is completely captured by the evolution 
(heat) kernel K(g,g',s) = K(gg_x,s) [II, EH, in the 
sense that given the initial condition ip(go,0), we have: 
i>(9,s) = J dg o K(g,g o ,s)ip{g o ,0). The various phys- 
ical propagators are obtained from the evolution ker- 
nel according to how the proper time is integrated out; 
the Hadamard function is obtained via the expression: 
H(g,g',m 2 ) = -i /+~ dsK(g,g> \ S ) e -* m2s , while re- 
stricting the range of integration to positive proper times 
only gives the Feynman propagator Gp(g,g' ,m 2 ) = 



ij °° dsK(g,g',s)e- 



prescription to 



m s , where the usual Feynman 
-ieisassumedforconvergence.lt 
turns out that the functions entering the expressions 
for the quantum amplitudes of all current spin foam mod- 
els correspond to the Hadamard 2-point functions for a 
scalar field on the relevant Lie group/homogeneous space 
with to 2 = —C(pf) where C(pf) is the Casimir eigenvalue 
of the representation labelling the face of the 2-complcx 
(therefore the link of the spin network whose evaluation 
gives the amplitude for the vertex): 



H(9 f ,m 2 ) = —^A Pfev D p ^(6 f ) 



+ OC 



= -i / ds nv K M (d Sl T v )e +lC{p ^ T ' , 



with the "oriented proper time" being T = p v s, A p is 
the dimension of the representation p in the Rieman- 
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nian case, or in the Lorentzian case (where the uni- 
tary representations are infinite dimensional) the contri- 
bution of the representation to the Plancherel measure, 
and the Casimir eigenvalues are C(p) = 2j(2j + n — 2) 
with j half-integer in the Riemannian case, and C{p) = 
+p 2 + (n — 2/2) 2 with p positive real in the Lorentzian 
case based on the timelike hyperboloid. Again a similar 
formula holds for the Ponzano-Regge models and for the 
models based on the (n-l)-dimensional DeSitter space. 
It is clear that the result is independent of the value of 
the various orientation data. On the other hand simply 
imposing the above restriction in the proper time integra- 
tion (that amount indeed to a causality restriction for the 
particle evolution Q), one obtains the expression for the 
Feynman propagator that one needs to define orientation 
dependent or causal spin foam models, and its expression 
in terms of representation functions of the 2nd kind, con- 
firms their interpretation as Wightman functions: 

r+ao 

G(tf/,m 2 ,^) = -i / ds flv K M (tif,Tf ]v )e +iC ^ T ^ = 



9(p v )E p i (d f ) + 8(-p v )E- Pf {d f ) 



= &V^VMw + @f,pf) + 0(-ih,)w-(#f,Pf)] ■ 

Similarly for the other spin foam models 0j- This is 
a non-trivial function of the orientation data, with the 
usual "time ordering" being replaced by a "timeless or- 
dering" ! 



V. A PROPER TIME FORMALISM AND THE 
FEYNMAN PROPAGATOR FOR SPIN FOAM 
QUANTUM GRAVITY 

Taking seriously this particle analogy and the associ- 
ated spin network feynmanology, and having also in mind 
the group field theory approach 0, we can generalise 
the current formulation of spin foam models to include 
a proper time variable, and obtaining a general expres- 
sion from which both orientation independent and causal 
models can be derived, simply changing the integration 
contour in proper time. The expression for the models 
dealt with above looks as follows: 

z =nn/ ds /i»nEnn / dx e A{p f ,x e ,T flv ) 

« f\v JC f PS v e\v JH - 

A(p f ,X e ,Tf\ v ) = H f Af(pf)U e Ae(fif\e) 

A similar formula holds for all the other models. This 
generalised expression encompasses both the usual un- 
oriented models and the new causal ones; indeed the first 
are obtained by choosing the extended range of integra- 
tion C — ( — oo, +oo) for the T variable (and this erases 
the dependence on p v and a e , while the quantum grav- 
ity Feynman propagator is obtained with C — (0, +oo). 



In this last case, a regularization prescription for conver- 
gence is implicit for each variable (so that the expression 
has to be understood in the complex domain): pf — > 
Pf +ie, i9y — > i)f +i5. The explicit form of the evolution 
kernel K differs of course in the different models and af- 
fects the exact form of the amplitudes [!j , that however 
all share the general structure here presented [25]. Let 
us discuss some properties of these new kind of transition 
amplitudes. First, all these models can be recast in the 
form of quantum causal histories models 0,^3; therefore 
they define quantum amplitudes for causal sets, if the un- 
derlying oriented graph does not contain closed timclike 
loops, and with the additional restriction of the vertices 
being (n + l)-valent in n dimensions. Second, while the 
un-oriented models can be related to the classical Regge 
action only in a asymptotic limit, when their vertex am- 
plitudes result in being proportional to the cosine of it 
0, HH Hflj l , here the connection with the Regge action is 
manifest without any approximation; for example, in the 
4-dimensional case the relevant evolution kernel has the 
form 14]: T) = 

causal vertex amplitude takes the form: 



e i 4T it an j j-^e 



A,. 



n 



-2mpf 



dsf\ v K{$f,Tf\ v )e 



iC{p t )Tj 



f P f smn $ / 

and therefore the product over vertex amplitudes in the 
spin foam model gives 



nA c =(nn^) 

V \ V f£v 11 1 J 

nn- 



I 



sinh 



i.e. the exponential of the Regge action (in first or- 
der formalism 0) for simplicial gravity, with an addi- 
tional contribution to the overall measure, as one would 



[25] The technology related to quantum particles on homogeneous 
spaces could have been used to generalise and then modify also 
the face and edge amplitudes. The reason why we have not done 
so is twofold: on the one hand the usual interpretation of these 
contributions to the spin foam models is that of a contribution 
to the overall measure, therefore the implementation of causality 
may be needed only at the vertex level, that is instead supposed 
to encode the dynamics of the theory; on the other hand, the 
form of the edge amplitudes in some version of the 4-d Barrett- 
Crane model is understood a s ar ising directly from the form of 
boundary spin network states and we prefer, at this stage, to 
keep that structure without modification. We believe, however, 
that alternative formulations of the models and possible defini- 
tions of modified spin network states, maybe in order to induce 
an orientatation-dependence in their structure, deserve further 
analysis. 
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expect from a sum-over-histories formulation of quan- 
tum gravity based on a simplicial discretization. The 
same is easily shown in the 4d Riemannian case and in 
the 3-dimensional case, while the proof for the higher 
dimensional models is less straightforward [9j- Third, 
the causal/oriented models seem to solve the problem, 
present in all BF-type formulations of quantum grav- 
ity [20I l2l| , of different isomorphic sectors related by a 
change of orientation all summed over in the path inte- 
gral quantization and thus interfering, leading to a dis- 
crepancy with the straightforward quantization of GR, 
and affecting the computation of geometric quantititcs 
in spin foam models |20f: in the causal models it seems 
instead that a restriction to the GR sector is achieved. 



VI. CONCLUSIONS 

Let us summarise our results: 1) we have linked the 
notion of causality in a Lorentzian context with the ori- 
entation of the 2-complcx on which spin foam models 
are based, and have identified the relevant data charac- 
terizing this orientation; 2) we have shown that all cur- 
rent spin foam models are orientation-independent; 3) 
using the technology of evolution kernels for quantum 
fields/particles on Lie groups/homogeneous spaces, we 
have constructed a generalised version of spin foam mod- 
els, introducing an extra variable with the interpretation 
of proper time; 4) we have proven that different ranges 
of integration for this proper time variable lead to differ- 
ent classes of spin foam models: one corresponds to the 
usual ones, interpreted as the quantum gravity analogue 
of the Hadamard function of QFT or, equivalenty, in a 
canonical interpretation, as a covariant definition of the 



inner product between quantum gravity states; the other 
is a new class of models and corresponds to the quantum 
gravity analogue of the Feynman propagator in QFT, i.e. 
a causal transition amplitude, non-trivial function of the 
orientation data, that implies a notion of "timeless or- 
dering" ; 5) we have shown how the causal model is man- 
ifestly related to simplicial gravity in the 4d Lorentzian 
case. All these results hold true in full generality, for 
the type of spin foam models considered, i.e. regardless 
of the spacetime dimension and signature. These results 
open quite a few lines of possible further research. We list 
some of them: 1) a group field theory formulation of gen- 
eralised and oriented models (work on this is in progress); 
2) an investigation of the possibility of defining a notion 
of "positive and negative energy" sectors, also at the level 
of spin network states, based on orientation/causal prop- 
erties, in a timeless framework; 3) the study of the role 
and significance of the proper time formalism, and of the 
possible link between this group theoretic proper time 
parameter and the lapse function of canonical quantum 
gravity; 4) an analysis of the relationship between the 
spin foam approach and other approaches to quantum 
gravity, starting from the new causal models, that seem 
to be the at the point of convergence of simplicial quan- 
tum gravity dynamical triangulations (when a sum 
over 2-complexes is implemented) |23| and causal sets, 
in addition to canonical loop quantum gravity |24j : 5) 
the new models can be interpreted as defining the ma- 
trix elements of an "evolution operator" , the basis for 
any "scattering calculation" , whose property could be 
studied, also in order to understand whether a notion 
of unitary evolution is feasible in Quantum Gravity, in 
absence of an external time coordinate, but with a clear 
notion of causality. 
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